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Abstract. We discuss the process of opinion formation in a completely 
homogeneous, democratic population using a class of probabilistic cellu- 
lar automata models with two absorbing states. Each individual can have 
one of two opinions that can change according to that of his neighbors. 
It is dominated by an overwhelming majority and can disagree against 
a marginal one. We present the phase diagram in the mean field ap- 
proximation and from numerical experiments for the simplest nontrivial 
case. For arbitrarily large neighborhoods we discuss the mean field re- 
sults for a non-conformist society, where each individual adheres to the 
overwhelming majority of its neighbors and choses an opposite opinion in 
other cases. Mean field results and preliminary lattice simulations with 
long-range connections among individuals show the presence of coherent 
temporal oscillations of the population. 



1 Modeling social pressure and political transitions 

What happens to a society when a large fraction of people switches from a con- 
formist to a non-conformist attitude? Is the transition smooth or revolutionary? 
These important questions, whose answers can make the difference between two 
well-known political points of view, is approached using a theoretical model, in 
the spirit of Latane's social impact theory [1,2]. 

We assume that one's own inclination towards political choices originates 
from a mixture of a certain degree of conformism and non-conformism. Con- 
formists tend to agree with the local community majority, that is with the av- 
erage opinion in a given neighborhood, while non-conformists do the opposite. 
However, an overwhelming majority in the neighborhood (which includes the 
subject itself) is always followed. 

We shall study here the case of a homogeneous population, i.e. a homoge- 
neous democratic society. 1 It may be considered as the annealed version of a 

1 The case with strong leader was studied in Ref [3]. 



real population, which is supposedly composed by a mixture of conformist and 
non-conformist people who do not change easily their attitude. 

In Sec. 2 we introduce a class of probabilistic cellular automata characterized 
by the size 2r+l of the neighborhood, a majority threshold q, a coupling constant 
J and an external field H . 

We are interested in the two extreme cases: people living on a one-dimensional 
lattice, interacting only with their nearest neighbors (r = 1) and people inter- 
acting with a mean-field opinion. 2 

In Sec. 3 we present the simplest case where each individual interacts with his 
two nearest neighbors (r = 1), the mean field phase diagram and the one found 
from numerical experiments. For this simple case, we find a complex behavior 
which includes first and second order phase transitions, a chaotic region and 
the presence of two universality classes [6] . In Sec. 4 we discuss the mean field 
behavior of the model for arbitrary neighborhoods and majority thresholds when 
the external field is zero and the coupling constant is negative (non-conformist 
society). The phase diagram of the model exhibits a large region of coherent 
temporal oscillations of the whole populations, either chaotic or regular. These 
oscillations are still present in the lattice version with a sufficient large fraction 
of long-range connections among individuals, due to the small- world effect [7]. 

2 The model 

We denote by x\ the opinion assumed by individual i at time t. We shall limit 
to two opinions, denoted by —1 and 1 as usual in spin models. The system is 
composed by L individuals arranged on a one-dimensional lattice. All operations 
on spatial indices are assumed to be modulo L (periodic boundary conditions). 
The time is considered discontinuous (corresponding, for instance, to election 
events). The state of the whole system at time t is given by a;' = (xq, . . . , x t L _ 1 ) 
with x\ G {-1, 1}; 

The individual opinion is formed according to a local community "pressure" 
and a global influence. In order to avoid a tie, we assume that the local commu- 
nity is formed by 2r+ 1 individuals, counting on equal ground the opinion of the 
individual himself at previous time. The average opinion of the local community 
around site i at time t is denoted by m* = X)j=-r x i+j- 

The control parameters are the probabilities p s of choosing opinion 1 at time 
t + 1 if this opinion is shared by s people in the local community, i.e. if the local 
"field" isTO = 2s-2r-l. 

Let J be a parameter controlling the influence of the local field in the opinion 
formation process and H be the external social pressure. The probability p s are 
given by 

Ps = P( m +2r+l)/2 °C exp(iJ + Jm). 

One could think to H as the television influence, and J as educational effects. 
H pushes towards one opinion or the other, and people educated towards con- 

2 Related models in one and two dimensions have been studied in Refs [4, 5]. 



formism will have J > 0, while non-conformists will have J < 0. In the statistical 
mechanics lingo, all parameters are rescaled to include the temperature. 

The hypothesis of alignment to overwhelming local majority is represented 
by a parameter q, indicating the critical size of local majority. If s < q (m < 
2q-2r- 1), then = -1, and if s > 2r + 1 - q (m > 2r + 1 - 2q), then 
** +1 = l- 

In summary, the local transition probabilities of the model are 

!0 if s < q; 
AB S 

TTAB* iS 9<s<2r + l-q; (1) 

1 if s > 2r + 1 - q; 

where A = exp[2ff + 2J(2r - 1)] and B = exp(4J). 

For q = the model reduces to an Ising spin system. For all q > we 
have two absorbing homogeneous states, x = — 1 (c = 0) and x = 1 (c = 1) 
corresponding to infinite coupling (or zero temperature) in the statistical me- 
chanical sense. With these assumptions, the model reduces to a one-dimensional, 
one-component, totalistic cellular automaton with two absorbing states. 

The order parameter is the fraction c of people sharing opinion 1. 3 It is 
zero or one in the two absorbing states, and assumes other values in the active 
phase. The model is symmetric since the two absorbing states have the same 
importance. 



3 A simple case 

Let us study the model for the simplest, nontrivial case: one dimensional lattice, 
r = 1 and q = 1. We have a probabilistic cellular automaton with three inputs, 
and two free control parameters, p\ and P2, while po = and p 3 = 1, according 
to Eq. (1). 

One can also invert the relation between (H, J) and {p\,pi): 



4 °v(i-pi)(i-P2)y 4 °\pi(i-p 2 ) 

The diagonal p 2 = p\ corresponds to J — and the diagonal pi = 1 — p\ to 
H = 0. 

At the boundaries of probability intervals, we have four deterministic (ele- 
mentary) cellular automata, that we denote T3, T23, T13 and T123, where the 
digits indicate the values of s for which p s = 1 [8] . 

Rule T3 (T123) brings the system into the c = (c = 1) absorbing state 
except for the special case of the initial configuration homogeneously composed 
by the opposite opinion. 

Rule T23 is a strict majority rule, whose evolution starting from a random 
configuration leads to the formation of frozen patches of zeros and ones in a 

3 The usual order parameter for magnetic system is the magnetization M = 2c — 1. 
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Fig. 1. Mean-field phase diagram for the density c coded as gray levels ranging from 
white (c = 0) to black (c = 1). The dashed upper-left region denotes the coexistence 
phase, in which both states c = and c = 1 are stable, and the final state depends on 
the initial density (first-order transition) . 



few time steps. A small variation of the probabilities induces fluctuations in the 
position of the patches. Since the patches disappear when the boundaries collide, 
the system is equivalent to a model of annihilating random walks, which, in one 
dimensions, evolves towards one of the two absorbing states according to the 
asymmetries in the probabilities or in the initial configuration. 

Rule Tf 3 on the other hand is a "chaotic" one, 4 leading to irregular pattern 
for almost all initial conditions (except for the absorbing states) . These patterns 
are microscopically very sensitive to perturbations, but very robust for what 
concerns global quantities (magnetization). 

The role of frustrations (the difficulty of choosing a stable opinion) is eviden- 
tiatcd by the following considerations. Rule T23 corresponds to a ferromagnetic 
Ising model at zero temperature, so an individual can simply align with the local 
majority with no frustrations. On the contrary, rule T13 is unable to converge 
towards an absorbing state (always present), because these states are unstable: 
a single individual disagreeing with the global opinion triggers a flip in all local 
community to which he/she belongs to. It is possible to quantify these concepts 
by defining stability parameters similar to Lyapunov exponents [10]. 

We start by studying the mean-field approximation for the generic case, with 
Pi and P2 different from zero or one. 

Let c and c' denote the density of opinion 1 at times t and t + 1 respectively. 
We have 

c' = 3pic (1 - cf + 3p 2 c 2 (1 - c) + c 3 . 




4 Called rule 150 in Ref. [9] 



Fig. 2. Numerical phase diagram for the density c (left) and hysteresis region for several 
values of the noise e and relaxation time T (right). Same color code as in Figure 1. 



This map has three fixed points, the state x = — 1 (c = 0), the state x = 1 
(c = 1) and a microscopically disordered state (0 < c < 1). The model is 
obviously symmetric under the changes pi — > 1 — p2, P2 — > 1 — Pi and x — > 1 — x, 
implying a fundamental equivalence of political opinions in this model. The 
stability of fixed points marks the different phases, as shown in Fig. 1. 

The stability of the state c = 1 (c = 0) corresponds to large social pressure 
towards opinion f (— 1). The value of J determines if a change in social pressure 
corresponds to a smooth or abrupt transition. 



3.1 Phase transitions and damage spreading in the lattice case 

The numerical phase diagram of the model starting from a random initial state 
with c° = 0.5 is shown in Fig. 2. The scenario is qualitatively the same as 
predicted by the mean- field analysis. In the upper- left part of the diagram both 
states c = and c = 1 are stable. In this region the final fate of the system 
depends on the initial configuration. 

Due to the symmetry of the model, the two second-order phase transition 
curves meet at a bicritical point (p t , 1 — p t ) where the first-order phase transition 
line ends. Crossing the second-order phase boundaries on a line parallel to the 
diagonal p\ = p2, the density c exhibits two critical transitions, as shown in the 
inset of the right panel of Fig. 2. Approaching the bicritical point the critical 
region becomes smaller, and corrections to scaling increase. Finally, at the bicrit- 
ical point, the two transitions coalesce into a single discontinuous (first-order) 
one. 

First-order phase transitions are usually associated to a hysteresis cycle due 
to the coexistence of two stable states. To visualize the hysteresis loop (inset of 
the right panel of Fig. 2) we modify the model slightly by letting po = 1 — ps = e 



with e < 1. In this way the configurations x = — 1 and x = 1 are no longer 
absorbing. This brings the model back into the class of equilibrium models for 
which there is no phase transition in one dimension but metastable states can 
nevertheless persist for long times. The width of the hysteresis cycle, shown in 
the right panel of Fig. 2, depends on the value of e and the relaxation time T. 

We study the asymptotic density as p\ and P2 move on a line with slope 
1 inside the dashed region of Fig. 1. For p\ close to zero, the model has only 
one stable state, close to the state c = 0. As pi increases adiabatically, the new 
asymptotic density will still assume this value even when the state c = 1 become 
stable. Eventually the first state become unstable, and the asymptotic density 
jumps to the stable fixed point close to the state c = 1. Going backwards on 
the same line, the asymptotic density will be close to one until that fixed point 
disappears and it will jump back to a small value close to zero. 

Although not very interesting from the point of view of opinion formation 
models, the problem of universality classes in the presence of absorbing states 
have attracted a lot of attention by the theoretical physics community in recent 
years [11, 12]. For completeness we report here the main results [6]. 

It is possible to show that on the symmetry line one can reformulate the 
problem in terms of the dynamics of kinks between patches of empty and oc- 
cupied sites. Since the kinks are created and annihilated in pairs, the dynamics 
conserves the initial number of kinks modulo two. In this way we can present an 
exact mapping between a model with symmetric absorbing states and one with 
parity conservation. 

Outside the symmetry line the system belongs to the directed percolation 
universality class [13]. We performed simulations starting either from one and 
two kinks. In both cases p t = 0.460(2), but the exponents were found to be 
different. Due to the conservation of the number of kinks modulo two, starting 
from a single site one cannot observe the relaxation to the absorbing state, and 
thus S = 0. In this case r\ — 0.292(5), z — 1.153(5). On the other hand, starting 
with two neighboring kinks, we find r] = 0.00(2), S = 0.285(5), and z = 1.18(2). 
These results are consistent with the parity conservation universality class [3,4]. 

Let us now turn to the sensitivity of the model to a variation in the initial 
configuration, i.e. to the study of damage spreading or, equivalently, to the lo- 
cation of the chaotic phase. Given two replicas x and y, we define the difference 
w as w = x (B y, where the symbol denotes the sum modulus two. 

The damage h is defined as the fraction of sites in which w = 1, i.e. as the 
Hamming distance between the configurations x and y. We study the case of 
maximal correlations by using just one random number per site, corresponding 
to the smallest possible chaotic region [14]. 

In Fig. 3 the region in which the damage spreads is shown near the lower- 
right corner (chaotic domain). Outside this region small spots appear near the 
phase boundaries, due to the divergence of the relaxation time (second-order 
transitions) or because a small difference in the initial configuration can bring 
the system to a different absorbing state (first-order transition). The chaotic 
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Fig. 3. Phase diagram for the damage found numerically by considering the evolution 
starting from uncorrelated configurations with initial density equal to 0.5. 



domain is stable regardless of the initial density. On the line P2 = the critical 
point of the density and that of the damage spreading coincide. 



3.2 Reconstruction of the potential 

An important point in the study of systems exhibiting absorbing states is the 
formulation of a coarse-grained description using a Langevin equation. It is gen- 
erally accepted that the directed percolation universal behavior is represented 

by ' 

dc(x, t) 



dt 



= ac(x, t) — bc 2 (x, t) + \7 2 c(x, t) + \J c{x, t)a(x, t), 



where c(x, t) is the density field, a and b are control parameters and a is a Gaus- 
sian noise with correlations (a(x, t)a(x',t')) = 5 x _ x <5 t ,t'- The diffusion coefficient 
has been absorbed into the parameters a and b and the time scale. 

It is possible to introduce a zero-dimensional approximation to the model by 
averaging over the time and the space, assuming that the system has entered a 
metastable state. In this approximation, the size of the original systems enters 
through the renormalized coefficients o, b, 

dc(x,t) _ 



dt 



= ac(x, t) — bc 2 (x, t) + \J c(x, tja(x, t), 



where also the time scale has been renormalized. 
The associated Fokkcr-Planck equation is 

dP(c,t) d . t 2\ t->/ \ 19 2 n , x 
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Fig. 4. Reconstruction of potential V(c) for P2 = (left) and for the kink dynamics 
on the line P2 — 1 — Pi (right). 

where P(c, t) is the probability of observing a density c at time t. One possible 
solution is a <5-peak centered at the origin, corresponding to the absorbing state. 

By considering only those trajectories that do not enter the absorbing state 
during the observation time, one can impose a detailed balance condition, whose 
effective agreement with the actual probability distribution has to be checked a 
posteriori. The effective potential V is defined as V(c) = — log(P(c)) and can be 
found from the actual simulations. 

In the left panel of Fig. 4 we show the profile of the reconstructed potential 
V for some values of p around the critical value on the line p>i = 0, over which 
the model belongs to the DP universality class. One can observe that the curve 
becomes broader in the vicinity of the critical point, in correspondence of the 
divergence of critical fluctuations x ~ \p ~~ Pel _7 > l' — 0-54 [15]. By repeating 
the same type of simulations for the kink dynamics (random initial condition) , 
we obtain slightly different curves, as shown in the right panel of Fig. 4. We 
notice that all curves have roughly the same width. Indeed, the exponent 7' for 
systems in the PC universality class is believed to be exactly [16], as given 
by the scaling relation 7' = dv± — 2/3 [15]. Clearly, much more information can 
be obtained from the knowledge of P(c), either by direct numerical simulations 
or dynamical mean field trough finite scale analysis, as shown for instance in 
Ref. [17]. 

4 Larger neighborhoods 

In order to study the effects of a larger neighborhood and different threshold 
values q, let us start with the well known two-dimensional "Vote" model. It is 
defined on a square lattice, with a Moore neighborhood composed by 9 neigh- 
bors, synthetically denoted M in the following [8]. If a strict majority rule q = 4 
is applied (rule M56789, same convention as in Sec. 2) to a random initial con- 
figuration, one observes the rapid quenching of small clusters of ones and zero, 
similar to what happens with rule T23 in the one-dimensional case. A small noise 
quickly leads the system to an absorbing state. On the other hand, a small frus- 
tration q = 3 (rule M46789) for an initial density c° = 0.5 leads to a metastable 
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Fig. 5. Case H — 0, J = -co: the mean field map Eq. (2) (left) for r = 10 and q — 2 and 
the mean-field r — q phase diagram (right). In the phase diagram the absorbing states 
are always present. Points mark parameter values for which the absorbing states are 
the only stable attractors. A plus sign denotes period-2 temporal oscillations between 
absorbing states, a star denotes the presence of a stable point at c = 0.5, a cross (circle) 
denotes period-two (four) oscillations between two non-zero and non-one densities, 
triangles denote chaotic oscillations. 



point formed by patches that evolve very slowly towards one of the two absorb- 
ing states. However, this metastable state is given by the perfect balance among 
the absorbing states. If one starts with a different initial "magnetization", one 
of the two absorbing phases quickly dominates, except for small imperfections 
that disappear when a small noise is added. 

In the general case, the mean-field equation is 

^=E( 2r , + 1 )^ 1 - c ) ar+1 "' p " (2) 

s=0 V S J 

sketched in the left panel of Fig. 5. 

We studied the asymptotic behavior of this map for different values of r and 
q. For a given r there is always a critical q c value of q for which the active phase 
disappears, with an approximate correspondence q c ~ 4/5r. The active phase is 
favored by the presence of frustrations and the absence of the external pressure, 
i.e. for J < and H = 0. We performed extensive computations for the extreme 
case J = — oo, H = corresponding to a society of non-conformists without 
television. As shown in the right panel of Fig. 5, by increasing the neighborhood 
size r, one starts observing oscillations in the equilibrium process. This is evident 
in the limit of infinite neighborhood: the parallel dynamics induced by elections 
(in our model) makes individual tend to disalign from the marginal majority, 
originating temporal oscillations that can have a chaotic character. Since the 
absorbing states are present, and they are stable for q > 0, the coherent oscil- 



lations of the population can bring the system into one absorbing state. This is 
reflected by the "hole" in the active phase in the mean field phase diagram. 

Preliminary lattice simulations (not shown) reveal that this behavior is still 
present if there is a sufficiently large fraction of long-range connections due to 
the small- world effect [7], while the active phase is compact in the short-range 
case. 

5 Conclusions 

Although this model is quite rough, there are aspects that present some analogies 
with the behavior of a real society. In particular, the role of education, repre- 
sented by the J parameter. In a society dominated by conformists, the transitions 
are difficult and catastrophic, while in the opposite case of non-conformist people 
the transitions are smooth. However, in the latter case a great social pressure is 
needed to gain the majority. 

On the other hand, if the neighborhood is large and non-conformism is 
present, one can observe the phenomenon of political instabilities given by tem- 
poral oscillations of population opinion, which could be interpreted as a symptom 
of healthy democracy. 
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